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Rotational quenching of H2CO by molecular hydrogen: 
cross-sections, rates, pressure broadening. 
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ABSTRACT 



Grenoble (IPAQ) UMR 527^, Grenoble F-38041, France 



We compute the rotational quenching rates of the first 81 rotational levels of 
ortho- and para-H2C0 in collision with ortho- and para-H2, for a temperature range 
of 10-300 K. We make use of the quantum close-coupling and coupled-states scattering 
methods combined with the high accuracy potential energy surface of Troscompt et 
al. (2009aL Rates are significantly different from the scaled rates of H2CO in collision 
with He; consequently, critical densities are noticeably lower. We compare a full close- 
coupling computation of pressure broadening cross sections with experimental data 
and sho w that our results are compatible with the low temperature measurements of 
Mengel & De Lucia (20001, for a spin temperature of H2 around 50 K. 
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1 INTRODUCTION 

In order to relate quantitatively the observed rotational 
spectra and the molecular abundances, knowledge of the 
relative importance of photon induced transitions and col- 
lision induced transitions is imperative. While the Einstein 
A and B coefficients can be obtained experimentally (see 
e.g. th e JPL and CDMS databases |Drouin| ( |2012[ ); |Miiller 
|et al.| ( |2005[ )), the rates of molecular coUisional excita- 
tion/quenching need to be calculated with help of precise 
microscopic frameworks. Many molecular quenching rates 
have been put forward in the last 40 years, either for molec- 
ular collisions with He or for collisions with H, electrons 
and H2. In many types of interstellar regions, especially so 
when molecular complexity is present, the main collider is 
molecular hydrogen, H2. Recently, a renewed large effort has 
been devoted to compute the rotational quenching rates of 
molecules by H2 (see references in the review^ by Van 
der Tak (20111) and most collisional data are avail- 
able in t he BASECOL [Dubernet et. aL\ ( |2012| ) and 



LAMDA Schoier et al. ( 2005 1 databases. Special em 



phasis was put on molecules specific to the Herschel Space 
Observatory, and among them, water. However, complex or- 
ganic molecules continue to play a prominent role in the un- 
derstanding of proto-stellar evolution as well as being probes 
for various interstellar environments. 

Among all organic molecules, formaldehyde (H2CO) is 
especially abundant, since it is the first stable molecule 
resulting from the hydrogenation of the ubiquitous CO 
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molecule ( Peters et al.|[20TT |. Being abundant and display- 
ing a large range of both transition frequencies and energy 
levels, formaldehyde is a tool of choice to probe the physical 



conditions of the gaseous interstellar matter (Maret et al. 
|2004[ [Kama et al.||2013"l ) . By using models that do not sup- 
pose Local Thermodynamical Equilibrium (LTE), it is pos- 
sible to reach reliable estimates of the molecular abundance 
of H2CO together with the other parameters of the gas, its 



temperature and density ( Ceccarelli et al. 2003 Van der Tak 



et al. 20071. The famous formaldehyde "anomalous" 



absorption was also shown to provide a probe of the 
ortho-to-para ratio of H2 (Troscompt et al. 2009b|) 
and a distance-independent tracer of the cosmic star 



formation history ( ,Darling &: Zeiger||20l2 |. 

Several computations of the excitation rates of H2CO 
have been proposed in the literature. Excitation by He 
atoms, being easier to perform, has been used for long as 
a model for excitation by H2, even if it is known for be- 



ing an underestimation of unknown precision ( Green et al 



1978 Green 1991 Sharma et al. 2012 1. In our group, we 



computed ( [Troscompt et al.||2009a[ ) the quenching rates of 
H2CO by H2, but only for the low levels of excitation of 
ortho-formaldehyde (only the first 10 levels) and for low ki- 
netic temperatures (T < 100 K). These computations were 
based on a high precision Potential Energy Surface (PES) 
for the van der Waals interaction H2CO — H2. In the present 
paper, we extend our previous computations to a much 
broader range of rotational energies and temperatures, for 
both ortho- and para-H2CO, using the same PES. Also, in 
order to assess the precision of this PES, we compare ex- 
perimental pressure broadening cross sections measured by 
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Mengel & De Lucia (20001 to our own computations, in a 



manner similar to our recent works on the rotational exci 
tation of H 2O and CO ( jDrouin fc Wiesenfeld|[2012b| [Faure 
et al.|2013| ). 

We organize the paper as follows. Section[2]describes the 
details of the scattering computation. Section [3] shows the 
results of our computations including cross sections, rates 
and pressure broadening. We end (section |4| with a discus- 
sion and a conclusion. 



2 SCATTERING CALCULATIONS 

As mentioned above, the PES for the interaction of 
H2CO and H2 has been described in our previous paper 
i Troscompt et al.|2009a l. The PES was computed for frozen 
monomer geometries (rigid-rotor approximation). The ge- 
ometries were those of the averaged distances and angles, 
at the ground vibrational state for both H2CO and H2. Us- 
ing average ground state geometries instead of equilibrium 
geometries have been shown to give more reliable results 
I Valiron et al.|2008 l. Since no modifications were per- 
formed on the PES, neither on the ab initio com- 
putations nor on the fits, the reader is deferred to 



Troscompt et al. ( 2009a I for all necessary details of the ah 
initio and fitting procedures (see also |Rist fc Faure| ( |2011[ ) 
for the fitting procedure). 

Both H2CO and H2 have each two identical ^H nuclei, 
of nuclear spin 1/2. Hence, both exist in para and ortho spin 
states. For H2, the para state (total nuclear spin 1 = 0) has 
even rotational states J2 = 0, 2, . . . (We denote by J2, the 
rotational quantum number of H2). The ortho state (/ = 1) 
on the opposite, has odd rotational states, J2 = 1,3,.... 
The H2 rotational constant is taken at _B = 60.853 cm~^. In 
the PES, the H-H distance is taken at 7? = 1.448761 A, its 



average value at ground vibrational level (Troscompt et al 



2009a 1. The H2CO molecule is an asymmetric rotor, with 
rotational constants (in cm"^) : A = 1.29543, B = 1.13419, 
C = 9.40552 and centrifugal parameters Dj = 0.251 10"^ 
DjK = 0.43110"", and Dk = 0.64810"^ . The rotational 
constant C is along the electric dipole, on the C2v symmetry 
axisPj Describing the rotational states by the usual Jk^Kq 
rotational (pseudo-)quantum numbers, para states (/ — 0) 
of H2CO correspond to Ka even and ortho states (I — 1) to 
Ka odd (the main quantum rotational number for H2CO is 
denoted by J throughout this paper). 

Scattering calculations were performed for all levels 
with a rotational energy Erot ^ 210 cm~^, that is up to 
JkaKc =1037 for ortho-H2CO and JkaKc = 744 for para- 
H2CO. We computed the inelastic cross sections with a col- 
lision energy 1 ^ -Ecoii ^ 1000 cm~^. This slightly extends 



the previous H2CO-He computations of Green (19911 and 



extends our previous H2CO-H2 computations by an order 
of magnitude. Experimental ( Bocquet et al.|1996 l and com- 
puted para and ortho rotational levels of H2CO are given 
in table 111 along with the level numbering used in the re- 
sults. Note that restricting energies to E i^ 210 cm~^ entail 
Ka ^ 4. Also, for high enough J values, pairs of succes- 
sive Ka Kc and Ka Kc + 1 levels tend to be degenerate. 



for all practical purposes. This degeneracy is all the more 
precise that the Kc values are small. Inspection of table [T] 
shows that the experimental - theoretical energy level dif- 
ferences AE are very small: \AE\ ^ 0.43 cm~^ (mean value 
{|A£;|) = 0.024 cm"^). 

All scattering calculations have been performed with 
the OpenMP version of the MOLSCAT cod(Q The reduced 
mass for H2CO-H2 is 1.889053 amu. The coupled-channel 
(CC) and coupled-states (CS) equations were integrated us- 
ing the diabatic modified log-derivative propagator. 

The rotational basis set is devised as follows. For all 
scattering energies, if Jk^Kc is the last open channels, 
the J' — J+l,J-l-2 values were added to the basis. 
However, since a given J rotational number spans a large 
amount of rotational energy, rotational levels were capped, 
at -Bmax = 400 cm~^ for collision energies Etot ^ 175 cm^^, 
and increasing progressively to -Bmax = 1200 cm^^ for i^tot ~ 
1000 cm" ^. These large values of i?max are needed 
to converge cross sections; a similar effect was pre- 
viously observed for methyl-formate colliding with 
Helium ( jFaure et ar]|2011[ ). 

The rotational basis for ortho-H2 is J2 = 1. It has 
been shown for many systems that including J2 = 3 in 
the basis does not have a noticeable influence for temper- 



atures as low as 300 K (e.g. Daniel, Dubernet fc Grosjean 



(2011 1). The rotational basis for para-H2 proved diffi- 
cult to settle. We were able to use the J2 = 0, 2 basis 
set for CC par a-H2 - ortho-H2CO colli sions, wth 
-Etot < 130 cm"^ ( [Troscompt et al.||2009a[ ). For para- 
H2CO, because of the level structure, the J2 =0,2 
basis for CC computations proved to be practi- 
cally impossible for Etot }t 50 cm"^. As a result, w^e 
resorted only to a J2 = basis, both for ortho- 
and para-H2CO, stretching the CC computations 
as high as possible and continuing with the CS ap- 
proximation. In order to assess the importance of 
the J2 = 2 channel, however, w^e complemented the 
ortho-H2CO data w^ith a full CS J2 = 0, 2 compu- 
tation up to Etot ~ 900 cm" ^, using a coarse energy 
grid. This allow^ed us to check i) that cross sections 
for transitions with J2 = 2 — > are negligible, li) that 
cross sections for transitions with J2 = 2 — > 2 are 
very similar to those with J2 = 1 — >^ 1 (as observed 
for other systems, see e.g. Daniel, Du bernet &: Gros-| 
jean (2011)) and iii) that the difference between the 
basis sets J2 = and J2 = 0, 2 decreases ■with increas- 
ing collision energy, from an average of ~ 30 % to 
below 10 %. 

Because of the large number of expansion terms of 
the potential function in the spherical harmonic basis 



(Troscompt et al. 2009a I, several strategies have been de- 



vised in order to converge the inelastic scattering computa- 
tions in a reasonable amount of time (We arbitrarily tried to 
limit ourselves to 72 hours of clock time, for one energy scat- 
tering point, on 12 CPU cores). The radial propagation used 
a step size parameter STEPS = 15 except at collision energies 
below 10 cm~^ where was STEPS progressively increased up 
to 50. Also, RMAX values were progressively increased in the 



^ This is the MOLSCAT convention, differing from usual spectro- 
scopic conventions by the ordering of the axes. 
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low-energy regime from default to 100. Other propagation 
parameters were taken as the MOLSCAT default values. 

For total energies above approximately 330 cm~^ the 
coupled channels (CC) approach, exact in the fully con- 
verged limit, proved to be impractical. We thus had to resort 
to the usual coupled states (CS) approximation, with all its 
shortcomings. Like in our recent HDO-H2 rate calculations 
I jWiesenfeld, Scribano fc Faure||2011[ ), we used an additive 
constant to scale appropriately the CS cross sections by their 
counterpart CC values. Overlap of the CC and CS calcula- 
tions showed the validity of these corrections. The statis- 
tics of difference S between the CC and CS cross- 
sections was also examined. We found the follow^ing 
values (in A , s, the standard deviation): {\5\) ^ 0.5, 
s{5) ^ 2. Also, for collisions with ortho-H2, the total angular 
momentum Jtot was not stepped by unit values, as is done 
usually. For CC calculation above -Bcoii > 200 cm~^, a step 
in Jtot value jstep — 4 was used, with a careful checking of 
the convergence of the jstep > 1 procedure. This procedure 
was used throughout the CS computations. 

The energy grid was chosen to guarantee a good de- 
scription of the resonances, including those pertaining to 
the high lying rotational states. Also, since we aim at rates 
for temperatures up to 300 K, a particular care was taken to 
ensure both some economy in the computational load and 
a good convergence of the rate computation. Let us recall 
that the quenching rate from state i to state j, kji-i{T) (in 
cm'^sec"^) is related to the inelastic cross section aji-i{E) 
(in A^, with E, the collision energy) by the well known Boltz- 
mann average : 



3 RESULTS AND DISCUSSION 

3.1 Cross sections 

The inelastic cross sections have a general shape that is sim- 
ilar to all earlier findings, for collisions of a molecule with 
H2. As usual, the inelastic scattering with J2 = may be 
markedly different from the scattering with J2 > 0. This was 
observed and thoroughly discussed for H2O and HDO scat- 



tering computations ( Dubernet et al 



2009 



fc Grosjean 20lT|[Wiesenfeld, Scribano fc Faure|2011||Faure 



Daniel, Dubernet 



et al. 20121, and observed for a wide range of other colli- 



sions, SO2 and CI atoms being recent examples (Cernicharo 
et al.|2011 Lique fc Alexander|2012 |. Experiments with H2 
colliding with water molecules also extensively confirm this 



difference prouin fc Wiesenfeld 2012b[ [Yang et al.||2011 |. 
The situation with formaldehyde colliding with H2, J2 = 
and J2 > is however less clear, as some ortho-H2 and para- 
H2 collisions are nearly identical, especially for small sections 
and large A J. Examples of astrophysical significative cross- 
sections o-(-Econ) are given in figure |4] Those three cases are 
representative of all these astrophysically relevant types of 
sections that we examined: o-(-Btot; J2 = 0) ^ a{Etot; J2 > 
0) even if the o-(_Etot; J2 = 0) displays a richer resonance 
structure, because of the disappearance of the supplemen- 
tary quantum coupling J -I- J2 = J12. In each case, we ob- 
serve that (y(E^o\\\ J2 = 1) ~ o-(£'coii; Ji = 2), in structure 
and magnitude. Remember, however, that the threshold for 
scattering with J2 = 2 is B(R2) x 4 = 243.4 cm"^ higher 
than for J2 = 1. A full and reliable computation of H2CO 
scattering with para-H2, J2 ^ 2 is thus impossible for all 
practical purposes, in the present computer configurations. 



kj^i(T) 



_8_J_ 



3.2 Rates 



Gj^,{E)Ee^p{-E/T)AE, (1) 



where T and E are expressed in the same units, and /i 
is the coUisional reduced mass. The probability density 
Eexp{—E/T) going very slowly down with energy, it is cus- 
tomary, for calculating kj^i{T = To), to compute crj^i{E) 
up to E r^ 10 To (e.g., [Dubernet eFall ( [200^ ). This ap- 
proach is prohibitively time-intensive for a heavy molecule 
like H2CO, even within the CS approximation. For high ly- 
ing rotational states, we inspected carefully the actual nu- 
merical convergence of equation IT] and stopped our energy 
grid as soon as the rate was saturated by 10% in general, 
and 25% for the highest levels, with E^p > 150 cm~^. 

The computation of the pressure broadening cross- 
sections ( [Wiesenfeld fc Faure||2010l [Drouin fc Wiesenfeld| 
2012b I necessitates a very fine energy grid and a good con- 



trol of the elastic cross sections, much more difficult to ob- 
tain than the corresponding inelastic sections. Hence, we had 
to re-calculate all the S'-matrices with collision energies from 
2 to 70cm~^.Then, we compute the a^^(T) in order to com- 
pare it with the experimental results of [Mengel fc De Lucia| 
(20001, in a temperature range of 10 — 30 K. Both collisions 
with para-H2 (basis set, J2 — 0,2) and ortho-H2 (basis set, 
J2 = 1) were performed on an identical fine energy grid and 
in the CC formalism. 



All quenching rates for all levels of tab le [T] a re computed 
for the same temperature grid as Green ( 1991 1, lOK ^ T ^ 



300 K. The full table is deposited in the LAMDA database 
(ISch oier et al.|2005[ ) and BASECOL database ( [Dubernet et7 
jal.||2012l, and may be asked to the authors. The rates with 



ortho-H2 are based on the J2 = 1 sections only. In all cases, 
rates with para-H2 are given as a Boltzman average over 
the populations of J2 = 0, 2, with the further approximation 
of a{Ecoii;J2 = 2) ~ cr(_Ecoii; J2 = 1), when necessary (see 
above) . The influence of the J2 = 2 initial states may indeed 
be very large, as was also observed, in another context, for 



the pressure broadening of II2O by II2 ( Drouin fc Wiesenfeld 
|2012b[ ). 

FigureHlcompares globally present critical densities and 
critical densities from iGregm (1991 1, for electric-dipole al- 
lowed transitions. We compare our present rates with para- 
H2 with the properly scaled rates with He. We have the 
following definition on the critical density n*: 



z^i'<i ^i'<- 



zZj<i kj^i{T) 



(2) 



where i is the level under scrutiny and j' denotes all levels 
connected by a radiative transition, while j spans all levels. 
Left panel shows the scatter of critical densities ratios at T = 
300 K, while on the right panel, all ratios are averaged and 
plotted against temperature. The scatter is moderate, with 
no ratio exceeding 3. The right panels shows the evolution of 
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Table 1. H2CO rotational levels, in cm ^. The experimental values are from |Bocquet et al.| | |l996| ; the MOLSCAT 
values are calculated via the rotational constants described in the text. 



Level J Ka 



Para H2CO 
Kc MOLSCAT 



Experimental Level J Kj\ 



Ortho H2CO 
Kq Molscat 



Experimental 



1 











0.00000 


0.00000 


1 


1 


1 


1 


10.53897 


10.53900 


2 


1 





1 


2.42961 


2.42960 


2 


1 


1 





10.70021 


10.70010 


3 


2 





2 


7.28640 


7.28640 


3 


2 


1 


2 


15.23672 


15.23690 


4 


3 





3 


14.56547 


14.56550 


4 


2 


1 


1 


15.72044 


15.72020 


5 


4 





4 


24.25953 


24.25970 


5 


3 


1 


3 


22.28172 


22.28220 


6 


5 





5 


36.35891 


36.35920 


6 


3 


1 


2 


23.24914 


23.24870 


7 


2 


2 


1 


40.04024 


40.04020 


7 


4 


1 


4 


31.67205 


31.67290 


8 


2 


2 





40.04262 


40.04260 


8 


4 


1 


3 


33.28429 


33.28350 


9 


3 


2 


2 


47.32780 


47.32780 


9 


5 


1 


5 


43.40518 


43.40670 


10 


3 


2 


1 


47.33970 


47.33970 


10 


5 


1 


4 


45.82316 


45.82200 


11 


6 





6 


50.85168 


50.85240 


11 


6 


1 


6 


57.47805 


57.48040 


12 


4 


2 


3 


57.04241 


57.04250 


12 


6 


1 


5 


60.86224 


60.86050 


13 


4 


2 


2 


57.07810 


57.07800 


13 


7 


1 


7 


73.88710 


73.89070 


14 


7 





7 


67.72387 


67.72520 


14 


7 


1 


6 


78.39720 


78.39490 


15 


5 


2 


4 


69.18225 


69.18240 


15 


3 


3 


1 


88.23817 


88.23820 


16 


5 


2 


3 


69.26541 


69.26520 


16 


3 


3 





88.23819 


88.23830 


17 


6 


2 


5 


83.74502 


83.74530 


17 


8 


1 


8 


92.62834 


92.63370 


18 


6 


2 


4 


83.91097 


83.91060 


18 


4 


3 


2 


97.95762 


97.95760 


19 


8 





8 


86.95982 


86.96210 


19 


4 


3 


1 


97.95777 


97.95780 


10 


7 


2 


6 


100.72797 


100.72860 


20 


8 


1 


7 


98.42283 


98.41970 


21 


7 


2 


5 


101.02562 


101.02500 


21 


5 


3 


3 


110.10856 


110.10860 


22 


9 





9 


108.54265 


108.54640 


22 


5 


3 


2 


110.10917 


110.10920 


23 


8 


2 


7 


120.12792 


120.12910 


23 


9 


1 


9 


113.69738 


113.70500 


24 


8 


2 


6 


120.62140 


120.62030 


24 


9 


1 


8 


120.93290 


120.92900 


25 


10 





10 


132.45497 


132.46091 


25 


6 


3 


4 


124.69196 


124.69200 


26 


9 


2 


8 


141.94120 


141.94321 


26 


6 


3 


3 


124.69380 


124.69380 


27 


9 


2 


7 


142.71099 


142.70931 


27 


10 


1 


10 


137.08956 


137.10020 


28 


4 


4 





155.16894 


155.16940 


28 


7 


3 


5 


141.70886 


141.70889 


29 


4 


4 


1 


155.16894 


155.16940 


29 


7 


3 


4 


141.71348 


141.71350 


30 


11 





11 


158.67956 


158.68851 


30 


10 


1 


9 


145.92017 


145.91530 


31 


10 


2 


9 


166.16374 


166.16690 


31 


8 


3 


6 


161.16034 


161.16051 


32 


10 


2 


8 


167.30709 


167.30440 


32 


8 


3 


5 


161.17049 


161.17059 


33 


5 


4 


1 


167.31421 


167.31461 


33 


11 


1 


11 


162.79999 


162.81450 


34 


5 


4 


2 


167.31421 


167.31461 


34 


11 


1 


10 


173.37619 


173.37041 


35 


6 


4 


2 


181.88960 


181.88989 


35 


9 


3 


7 


183.04741 


183.04781 


36 


6 


4 


3 


181.88961 


181.88989 


36 


9 


3 


6 


183.06768 


183.06790 


37 


12 





12 


187.20017 


187.21330 


37 


12 


1 


12 


190.82361 


190.84309 


38 


11 


2 


10 


192.79099 


192.79581 


38 


12 


1 


11 


203.29127 


203.28461 


39 


11 


2 


9 


194.42165 


194.41750 


39 


10 


3 


8 


207.37096 


207.37180 


40 


7 


4 


3 


198.89572 


198.89600 


40 


10 


3 


7 


207.40854 


207.40891 


41 


7 


4 


4 


198.89574 


198.89600 















the average (n*) as a function of T. For higher temperatures, 
the effect of H2 being different of He diminishes, since higher 
colhsions energies are more sensitive to the hard walls of the 
target, where the potential grows exponentially (repulsion of 
the wave functions). H2 and He become more similar, and 
the rates tend one towards the other. 



3.3 Pressure broadening 

In order to assess the reliability of the PES, we found it use- 
ful to compare measure d and computed pressure-broadening 



very limited set of H2CO pressure broadening data exists 



cross sections, a^^(T) ( Wiesenfeld fc Faure|2010 Drouin & 



Wiesenfeld|2012b Faure et al.|2013 1. The main advantage of 
pressure broadening is that experiments and computations 
held both absolute quantities, with no scaling involved, rend- 
ing the comparison very meaningful. Unfortunately, only a 



for an H2 buffer gas ( Mengel fc De Lucia_2000J , for very low 
temperatures. The experimental results for the 2i2 — )■ 813 
transition are depicted in figure [3l along with full CC cal- 
culations of cr^^(T). We see that the para-H2 and ortho-H2 
computation bracket the experimental values, which are very 
well simulated by an ortho- to para-H2 ratio (opr) corre- 
sponding to a pseudo-equilibrium at ~ 50 K, correspond- 
ing to a OPR of 0.27. That the OPR of H2 may vary during 
the coUisional cooling experiment has been proved with cell 



walls covered with amorphous water ( [Drouin fc Wiesenfeld 
2012b I. Nothing is known for formaldehyde, and discussion 
with the authors of [Mengel fc De Lucia| ( [2000[ ) could not 
settle the case. We remain thus with a good plausibility ar- 
gument, as long as the pressure broadening experiments with 
H2 are not fully characterized. 
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4 DISCUSSION - CONCLUSION 

It is important to know which errors are to be expected, and 
to have some clues on how these errors might influence as- 
trophysical modeling. We expect the error on the rates k{T) 
to be uniformly increasing from low lying levels to higher 
lying ones, and also from low temperature to high temper- 
ature. Quantifying this error is very risky, as errors may 
arise from all the phases of the rate computation: ab initio 
computation, fits and long distance behavior of the PES, 
convergence of the CC/CS procedures, convergence of the 
averaging procedure of the sections with collision energy. 

Internal consistency with our earlier approach 
(Troscompt et al. 2009a I show differences less than 



10 %, for J ^ 5, T ^ 50 K. This shows that the convergence 
error in this domain cannot exceed 20-30%. In the high 
temperature domain (T ^ 150 K), it is safe to 
assume errors much larger than 30%, because of the 
poor convergence of formula nt, but still within 
a factor 2. The magnitude of the 5 values as well 
as the similarities betw^een the J2 = 1 — >^ 1 and 
J2 — 2 ^ 2 plead in favor of such a conservative 
value, as does the convergence of critical densities 
between our w^ork and the previous collisions with 
He by [Green] ([1991). 

Accuracy of the PES is very difficult to assess, without 
any firm experimental comparison, like has been done for 
H2O and to a lesser extent CO (Yang et al.|2011 Drouin & 



|Wiesenfeld|2012b[[Faure et al.|26l3[|Chefdeville et al.|2012 
Our results are, however, compatible with measurements of 



Mengel & De Lucia ( 2000 1. Because of the great importance 



and ubiquity of the formaldehyde molecule, further exper- 
iments would give indications on the precision of the PES 
and convergence procedures used in this paper. 

The relevance of the ortho-to-para ratio of H2 has been 
stressed several times already, as it may be of crucial impor- 
tance in order to correctly model the astrophysical environ- 
ments. While the difference in rates for ortho-H2 and para- 
H2 is large for low lying levels at low energies, see |Troscompt| 
et al. ( 2009b I for an application, this difference decreases at 



larger transition energies, as the hard walls of the PES play 
a more important role than the long range behavior. Indeed, 
the main increase in a{J2 > 0) with respect to (j(J2 = 0) 
is due to the averaging out of the quadrupolar moment and 
dipolar polarizability of II2 in its ground rotational states. 
The same is true, up to a global scale, for the difference in 
behavior between He and ortho-H2. 

Extension of these computations to higher levels and 
higher temperatures is by no means a difficult task, on the 
physics point of view, because of the rigid rotor structure of 
H2CO, its first bendin g frequency (the out-of-pl ane bend) 
arising at 1167.3 cm~^ ( Clouthier & Ramsay 1983 1. The true 



limiting factor arises from numerical load, with very large 
N X N matrices to invert and propagate (TV > 2000). The 
same situation arises for the excitation of heavier complex 
organic molecules, like methyl formate or dimethyl ether, 
which display many spectral lines very far from LTE, in 
various spectral surveys like Caux et al. (20111. Unfortu- 



nately, quasi classical trajectories methods are limited for 
all those cases because of ambiguities arising in the subse- 



quent quantization of rotational levels ( Faure & Wiesenfeld 



2004). Use of very large grid of computers and combinations 



of OpenMP/MPI approaches may overcome these difficul- 
ties. 

We have calculated an extensive set of low to medium 
temperature quenching rates, for all levels of H2CO below 
210 cm~^. These rates are ready to be incorporated in the 
various non-LTE models for the interstellar medium. At low 
temperatures, differences with earlier rates of H2CO collid- 
ing with Helium are very important and they remain notice- 
able at all temperatures, with ratios up to 50% at 300 K, 
where they are the most similar. Since [Troscompt et al.| 
(2009a) was tailored to be very precise at low tem- 
peratures for ortho-H2CO, we still recommend to 
use those rates for applications at T ^ 30 K. The 
present rates should have a large importance on the 
inferring of H2CO column densities, away from LTE 
conditions. 
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Figure 1. (Color online) Cross-sections (upper three panels) and rates (lower three panels) for three widely observed 



quenching transitions of ortho-H2CO (first column : 2ii 



li 



at Eu 



22.62 K; second column : 432 — > 831, at Eu 



140.93 K; third column : 735 —^ 634, at Eup = 203.9 K). In all panels the red (dashed) line denotes the H2 rotational J2 = 
initial and final state, the blue (solid) line, the same with J2 = 1 and the black (dashed-dotted), the same with J2 = 2. 
The red dash-dotted rates include the influence of the population of the J2 = 2 para state of H2 (see text). Rates and 
sections with J2 = •<-> J2 = 2 are two or three orders of magnitude lower and not depicted. Rates for collisions with Helium ( Gr een| 
|1991| l are the lowest rates in all three inferior panels (green color online). Note that for cross sections, for sake of clarity, the energy is 
the total energy. Note also that cross sections display usual resonance patterns which are not conspicuous at this log scale. 
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Figure 2. Left panel: ratio of critical densities (eq. [2|, present work to |Green| l |1991[ l properly scaled, for all levels, at T = 300 K. Right 
panel, average over all critical densities ratios, as a function of temperature. Both panels : open squares, ortho-H2CO with para-H2 
(J2 = 0) or He; filled diamonds para-H2CO with para-H2 (J2 = 0) or He. 
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Figure 3. (Color online) Pressure broadening of the 2i2 — > 813 
transition. Red symbols, measurements of |Mengel fc De Lucia] 
l |2000[ |. Dashed and dot-dashed lines, present computations with 
respectively ortho-H2 and para-H2. Full line, ortho-to-para of H2 
equivalent to a spin temperature of 50 K. 
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